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For the small-angle Bhabha-scattering process, we consider the error budget for 
the calculation of the LEP/SLC luminosity in the Monte Carlo event generator 
BHLUMI 4.04, from the standpoint of new calculations of exact results for the 
respective 0(a 2 ) photonic corrections in the context of the Yennie-Frautchi-Suura 
exponentiation. We find that an over-all precision tag for the currently available 
program BHLUMI 4.04 can be reduced from 0.11% to 0.061% at LEP1 and from 
0.25% to 0.122% at LEP2. For the large-angle Bhabha process, we present the 
Monte Carlo program BHWIDE and compare its predictions with predictions of 
other Monte Carlo programs as well as semi-analytical calculations. 



1 Introduction 

At LEP, for practical purposes, the process of Bhabha scattering, e + e~ — » 
e+e~, is divided into two classes depending on the kinematical regions^} the 
small-angle Bhabha (SABH) scattering, for the scattering angle 1° < e < 6°, 
and the large-angle Bhabha (LABH) scattering, for e > 10°. In general, this 
process is mediated by both the 7 and Z bosons, exchanged in both s- and 
t-channels. Thus, at the lowest order, there are 4 pure s- and ^-channel con- 
tributions and 6 interference terms (between 7 and Z and between the s- and 

a In the following we shall always consider the case with both the e+ and e~ detected. 
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^-channels) . At low angles, however, the Bhabha scattering is almost com- 
pletely dominated (> 99% of a cross section) by the pure QED process of the 
t-channel 7-exchange for which a very high accuracy in theoretical predictions 
can be achieved. For this reason SABH was chosen at LEP1 and LEP2 for the 
luminosity monitoring. At large angles, physical features of the Bhabha process 
at LEP1 and LEP2 are very different, as different Feynman-diagram contribu- 
tions dominate at these two energy regimes. At LEP1 energies, a dominating 
role is played by the s-channel Z-exchange, so this process is used, in parallel 
to other fermion-pair production, to determine properties of thc^Z-boson as 
well as to measure other important electroweak (EW) parameters^. At LEP2, 
LABH is dominated by the i-channel 7 exchange, so at these energies it resem- 
bles more SABH than the other two-fermion processes. Therefore, LABH is 
not very useful at LEP2 for testing the EW sector of the Standard Model (SM). 
It is considered mainly as a background for those processes that are of the main 
experimental interest at LEP2. The e + e _ -channel is investigated particularly 
in searches for possible "new physics" , such as SUSY, contact interactions, etc. 

Our discussion is organized as follows. In Section 2, we concentrate on 
the SABH scattering and discuss the new error budget for the Monte Carlo 
(MC) program BHLUMI 4.040 based on the exact 0{a 2 ) calculations. The 
LABH process is discussed in Section 3, where we give a brief description of the 
MC event generator BHWIDE0 and discuss its cross-checks and comparisons 
with those of other programs. Finally, Section 4 contains our conclusions and 
outlook. 

2 Small-Angle Bhabha Scattering 

Currently, new luminometers at LEP have given results on the luminosity 
process e + e~ — > e + e~ + n("f) at experimental precision tags below 0.1%. This 
should be compared with the prediction by the Krakow-Knoxville Collabo- 
ration in the Monte Carlo program BHLUMI 4.040, wherein the theoretical 
precision tag of 0.11% is realized for this process in the ALEPH SICAL-typeEl 
acceptance. If one combines the experimental results, one arrives at an ex- 
perimental precision of < 0.05%. Evidently, for the final EW precision tests 
data analysis for LEP1, it would be desirable to reduce the theoretical preci- 
sion tag on the luminosity cross section prediction, at least to the comparable 
0.05%- regime, in order not to obscure unnecessarily the comparison between 
experiment and the respective Standard Model of the electroweak interaction. 
With this aS|-pur primary motivation, we examined the error budget arrived 
at in Refs. 00 in view of recent exact results impacting both the technical 
and physical precision of the errors quoted in that budget. More precisely, if 
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one looks into the error budget shown in Table 1 of Rcf. □, one sees that the 
largest contribution is associated with the 0(a 2 ) photonic corrections, which 
contribute 0.1% in quadrature to the total 0.11% quoted for the total precision 
of the BHLUMI 4.04 prediction in these references for the ALEPHJjilCAL-type 
acceptance. Accordingly, we have used the exact results of Refs.&tM and the 
exact result of Ref . Ell to make a-more realistic estimate of the true size of this 
dominant error quoted in Refs.BQ. 

In re-examining the photonic corrections used in BHLUMI 4.04 at the 
0(a 2 ), which is the relevant order of the corrections, one needs to look at the 
approximations made in the matrix element used in the calculation encoded 
in the program in comparison to available exact results. This will allow us to 
re-assess the physical precision of the corresponding part of the BHLUMI 4.04 
matrix element, which is the exact 0(a 2 ) LL (leading- log) Yennie-Frautschi- 
Suura (YFS) exponentiated matrix clement. The phase-space integration of 
two hard photon emission in BHLUMI is exact (the LL approximations are 
only in the matrix element). Nevertheless, this four-body phase-space inte- 
gration should be cross-checked with another, independent, exact integration 
method. This check, which we have recently completed, will allow us to give 
a more realistic estimate of the technical precision of the realization of the 
corresponding aspect of the matrix element in BHLUMI 4.04. The net result 
is a new estimate of the total precision of the prediction of the luminosity cross 
section by BHLUMI 4.04 at LEP1 and LEP2 energies. 

Considering now the exact 0(a) correction to the single hard bremsstrah- 
lung in the luminosity process, we have implemented the results of Ref. into 
BHLUMI 4. xx and made a systematic study of the net change in the prediction 
for the luminosity relative to the prediction of BHLUMI 4.04 in which this 
correction is treated to the LL level. What we find is illustrated in Fig. [l]a 
for the ALEPH SICAL-type acceptance at the Z peak. In the language of 

~(2) 

the YFS theory, this correction enters the hard-photon residuals as p\ , the 
0{a 2 ) contribution to the one-hard-photon residual /3i. In Fig. |l|a, we show 
the difference between the corresponding LL result in BHLUMI 4.04 and: (1) 
our exact result as given in Ref. 0, (2) the approximate ansatz in Eq. (3.25) 
of Ref. 0, (3) the result (NLLB) of Ref. M, which is supposed to include the 
dominant non-LL effects, in ratio to the respective Born cross section. What 
we see is that the BHLUMI 4.04 results are within 0.02% of the exact result 
throughout the experimentally interesting regime 0.2 < 1 — z m i n < 1.0. This 
is the main reason why we will be able to reduce thfi estimated precision of the 
BHLUMI 4.04 prediction in comparison to Refs.00. 

Turning next to the technical precision of the 27 bremsstrahlung calcu- 
lation in BHLUMI 4.04, we have constructed a completely independent real- 
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Figure 1: (a) The pure second order Monte Carlo result for (3^ — P^Bhlumi differences 
for the SICAL Wide-Narrow trigger, divided by the Narrow-Narrow Born cross section; 
z m i n is as it is defined in Fig. 2 of Phys. Lett. B353 (1995) 362. (b) Comparison 
of Monte Carlo results for for the LL and exact matrix elements. The results are 
shown for the BHLUMI generator and for an alternative 'Test' generator for a technical 
precision test. 



ization of the two-photon phase-space integration compared to what is used 
in BHLUMI 4.04 by way of an independent MC algorithm. We have im- 
plemented this new MC realization of the exact two-photon phase-space and 
compared its result with that of BHLUMI 4.04's for the hard-photon residual 
$2 contribution to the luminosity cross section, both for the LL matrix ele- 
ment in BHLUMI 4.04 and for the exact matrix element for the two-photon 
bremsstrahlung of Refs.tffl. What we find is shown in Fig. for the ALEPH 
SICAL-type acceptance at the Z-peak. We find that the difference between 
the two realizations of the exact 27 bremsstrahlung is below 0.003% of the 
Born cross section. Moreover, we get an estimate of the physical precision of 
the LL approximation for this part of the cross section from compariruzJjhe LL 
and exact results as 0.012%, in agreement with our estimate in Refs.&H. 

Finally we turn to the exact result for the two-loop contribution of the 
hard-photon residual flo to the cross section in comparison to the LL result used 
for it in BHLUMI 4.04. We have analytically continued the result of Ref. E3 
from the s-channel to the t-channel for the required two-loop contribution to 
the respective charge form factor in QED. In this way using the YFS theory 
we have found that the difference between the LL result in BHLUMI 4.04 and 
the exact result corresponds to the shift of the function v in Eq. (2) of Ref. Il3 

by 

Aw « = (a ) 2 [(6 + 6C(3) - f - 4) L + 6 - 9«3) + (f - 2 In 2) ^ - , 
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where the big logarithm is defined as L = In \t\/m 2 and £(3) is the Riemann 
(■-function of the argument 3. For the ALEPH SICAL-type acceptance at the 
Z-peak, this corresponds to 0.014% in the cross section. 

Collecting the above results in quadrature, we obtain the result that the 
current calculation of the 0(a 2 ) photonic corrections in BHLUMI 4.04 are 
accurate to 0.027%. Using this result in Table 1 of Ref. □ we arrive at the 
precision tag 0.061% for the currently available calculation in BHLUMI 4.04 
at the Z peak. At the LEP2 energy of 176 GeV, if we repeat the analysis 
just described, we find that the corresponding precision of BHLUMI 4.04, for 
both the SICAL- and LCAL-type acceptances, inflow reduced to 0.122%, to 
be compared to the estimate of 0.25% in Refs. aB. The current situation is 
now summarized in Table [II A more detailed exposition of the results in this 
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Total 


0.11% 


0.061% 


0.25% 


0.122% 



Table 1: Summary of the total (physical+technical) theoretical uncertainty for a typical calori- 
metric detector. For LEP1, the above estimate is valid for the angular range within 1°— 3°, and 
for LEP2 it covers energies up to 176 GeV, and the angular ranges within l°-3° and 3°-6°. 



paper will appear elsewheret_3. 

Our result on the size of the error associated with the missing sub-leading 
bremsstrahlung correction at 0(a 2 ) in BHLUMI 4.04_which is 0.027%, agrees 
with the estimate of 0.03% made by Montagna et al.tB, using a structure func- 
tion convolution of a hard collinear external photon with an acollinear internal 
photon. As these authors have argued, while such a pairing of convolutions 
does not represent a complete set of photonic 0(a 2 L) corrections, one expects 
it to contain the bulk of such corrections. Indeed, our exact result of 0.027% 
shows that the approximation made in Ref.E3 does give the bulk of the respec- 
tive 0(a 2 L) correction. Evidently, the fact that we now have two independent 
results, one .exact, that presented by us in this paper, and one approximate, 
that in Ref. 113, which agree on the size of the error associated with the miss- 
ing photonic 0(a 2 L) correction in BHLUMI 4.04, enhances the results in this 
paper. 
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3 Large-Angle Bhabha Scattering 



At LEP1, the main physical quantities measured in LABH are: the total cross 
section er e and the forward-backward charge asymmetry Afb- A value of <j e is 
used to extract the partial Z decay width r e , while Afb is sensitive to the im- 
portant EW parameters, such as the top and Higgs masses. The experimental 
precision for LABH, after the final LEP1 data analysis is completed, is expected 
to be < 0.5% at the Z peak and ~ 1% at ±2 GeV off the peakH At LEP2, the 
e + e~-channel is considered mainly for the "new physical' searches, and there, 
the experimental precision is expected to reach 0.5%. Ell. On the theory side, 
several programs (bothlhe MC and semi-analytical ones) for LABH have been 
developed, see e.g. Ref.El, but a comprehensive analysis of the theoretical error 
is still missing^]. An important step in this direction was made during the '95 
Workshop "Physics at LEP2" H, where comparisons of several codes for LABH 
at LEP1 and LEP2 energies were performed. They showed, however, that the 
predictions of various programs can differ by as much as 2% at LEP1 and 4% 
at LEP2. In this section, we briefly describe our MC event generator for LABH 
called BHWIDE and discuss some important cross-checks of the program as 
well as comparisons of its predictions with the results of other programs. _ 
BHWIDE is based on the YFS exclusive exponentiation procedure E3, 
where all the IR singularities are summed-up to infinite order and cancelled 
out properly in the so-called YFS form factor. The remaining non-IR residuals, 
(3n \ corresponding to the emission of n-real photons, are calculated perturba- 
tively up to a given order I, where I > n, and (I — n) is a number of loops in 
the j3$ calculation. In BHWIDE an arbitrary number n of real photons with 
non-zero pr are generated according to the YFS MC method of Ref. a. The 
non-IR residuals (3$ are calculated up to 0(a), i.e. /3q and /3j correspond- 
ing to zero-real (one-loop) and one-real (zero- loop) photons, respectively, are 
included. In /3q we implemented—two libraries of the 0(a) virtual EW cor- 
rections: (1) the older one of Ref.tJ, which is not up, to date but can be useful 
for some tests, and (2) the more recent one of Ref. c3. When the genuine weak 
corrections are switched off (or numerically negligible) they are equivalent. In 
/?q we implemented two independent matrix elements for single-hard-photon 
radiation: (1) our calculation in terms of hclicity amplitudes, and (2) the 
formula of CALKULcj for the squared matrix element. We have checked that 
the above two representations agree numerically up to at least 6 digits on an 
event-by-event basis. This constitutes a very important technical cross-check 

''Recently, the analysis of the theoipfcical errors of two semi-analytical programs, ALIBABA 
and TOPAZ0, has been presented cj, but only for the LEP1 energies and for the so-called 
BARE acceptance. 
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of the implementation of the hard-photon matrix element in BHWIDE. 

The MC algorithm of BHWIDE is based on the algorithm of the program 
BHLUMI for SABHB with a few important modifications: (1) QED interfer- 
ences between the electron and positron lines ( "up-down" interferences) had to 
be reintroduced as they are important in LABH; (2) the full YFS form factor 
for the 2 — ► 2 process, including all s-, t- and it-channels, was implemented 13; 
(3) the exact 0(a) matrix element for the full BHABHA process was included. 
The multiphoton radiation is generated at the low-level MC stage as for the 
i-channel process, while the s-channel as well as all interferences are reintro- 
duced through appropriate MC weights. This means that the program is more 
efficient when the i-channel contribution is dominant, as e.g. at LEP2 energies; 
however, it proved to work well also at the Z peak. 

Having all necessary ingredients in the program and an appropriate algo- 
rithm for the MC event generation, we had to subject BHWIDE to several tests 
in order to check that it gives correct predictions. First, we wanted to know 
whether it reproduces the small-ansdc limit correctly. To this end we com- 
pared BHWIDE with BHLUMI 4.040 whose precision in the SABH regime, as 
shown in the previous section, is under the per-mille level. We found that for 
the angular acceptance of 1° < e < 10° the two programs agree within 0.1% 
(statistical error) for both the pure 0(a) QED corrections and the full YFS 
exponentiated cross sections at the energies of 5, 10, and 91.19 GeV. Then, 
we turned to large angles (40° < 9 e < 140°) and compared the pure QED 
(Z-exchange switched off) O(a) predictions of BHWIDE with the ones of the 
MC program OLDBISEj (a modernized version of the program X)LDBABlZI) 
whose technical precision was shown to be at the level of 0.02%Ej. We found 
an agreement between these two programs up to 0.05% (stat. error) for both 
the BARE and CALO acceptances as defined in Ref. El at the Z-peak energy. 
This should also remain valid at LEP2 energies, since without the Z contribu- 
tion there is no qualitative difference in LABH at these two energy regimes. 
The above result is a very important technical cross-check of BHWIDE, both 
for the implementation of the 0(a) QED corrections and for the correctness of 
the MC algorithms at large angles (the 0(a) result was extracted in BHWIDE 
from the full multiphoton YFS calculation). 

Finally, with all contributions/corrections at 0(a)YJ^ included, we com- 
pared BHWIDE with several MC and semi-analytical programs for LABH. 
These comparisons were first presented in Ref. El, and then updated in Ref. a. 
At LEP1 energies, the program TOPAZ0c3 is considered a semi- analytical 
benchmark for LABH. Its theoretical errors for the BARE-type acceptance, as 
recently estimated K are 0.4% (0.2%) at the Z peak and 0.3% (0.7%) on the 
wings (±2 GeV off the peak) for the maximum acollinearity of the final e + and 
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e~: 9™ a f = 10° (25°). From the tables and figures of RefsBB, we observe that 
for the BARE trigger BHWIDE agrees with TOPAZ0 within 0.4% (0.55%) at 
the Z peak and within 0.4% (0.95%) on the wings, for 9™f = 10° (25°). For 
the CALO trigger the agreement between these two programs is within 0.15% 
(0.25%) at the Z peak and 0.2% (0.55%) on the wings. We can see from the 
above results that the agreement between the two calculations improves con- 
siderably, as it should, for the CALO-type acceptance, which is closer to the 
real experiment. 

For LEP2 energies, the comparisons in Refs.0H were done for the CALO 
trigger only. Here, the discrepancies between various calculations are larger 
than at LEP1. This can be explained by the fact that most of the programs 
were constructed for the Z peak region, i.e. assuming that the s-channel 
contribution is dominant, while at the LEP2 energies this is not_the case. 
Here, a program designed for the i-chaimel process, as e.g. SABSPVE3, should 
be more reliable. The results in Refs.lM show that BHWIDE is within 1.5% of 
SABSPV for the whole LEP2 energy range. We have also checked that for the. 
BARE trigger BHWIDE agrees with the semi-analytical code ALIBABA c3 
within 0.3% at the pure 0(a) level and within 1% when all corrections are 
included, for the same energy range. 

From the above comparisons we can see that BHWIDE is in a better 
agreement with the semi-analytical benchmarks for LEP1 and LEP2 energies 
than any other MC event generator for LABH. A more detailed analysis of the 
BHWIDE theoretical precision is in progress now. 



4 Conclusions and Outlook 

We have discussed some aspects of the Bhabha scattering at LEP1 and LEP2. 
For the small-angle Bhabha process, we have presented the new error budget 
of the program BHLUMI 4.04 based on the exact 0(a 2 ) calculations. Wc 
have shown that the theoretical error for the luminosity measurement can be 
reduced now from 0.11% to 0.061% at LEP1 and from 0.25% to 0.122% at the 
LEP2 energy of 176 GeV. The predictions of BHLUMI 4.04 remain unchanged. 
The exact calculations can be included in the future version of the program if 
necessary. For the large-angle Bhabha process, we have presented the MC event 
generator BHWIDE and discussed some of its cross-checks and comparisons 
with other programs. From this we conclude that BHWIDE is the most precise 
MC event generator for LABH at LEP1 and LEP2. A comprehensive analysis 
of its theoretical errors is in progress. 
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